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Abstract —The uplink of interference-limited cellular networks 
with base stations that have large numbers of antennas and 
use linear Minimum-Mean-Square Error (MMSE) processing 
with power control is analyzed. Simple approximations, which 
are exact in an asymptotic sense, are provided for the spectral 
efficiencies (b/s/Hz) of links in these systems. It is also found 
that when the number of base-station antennas is moderately 
large, and the number of mobiles in the entire network is large, 
correlations between the transmit powers of mobiles within a 
given cell do not significantly influence the spectral efficiency of 
the system. As a result, mobiles can perform simple power control 
(e.g. fractional power control) that does not depend on other users 
in the network, reducing system complexity and Improving the 
analytical tractablllty of such systems. 

Index Terms —Massive MIMO, MMSE, power control. 

L Introduction 

Massive Multiple-Input Multiple-Output (MIMO) has been 
proposed as a technology to meet the exponential growth in 
demand for wireless communications that is forecasted for the 
near future. These systems use large antenna arrays at base 
stations to separate signals to/from users spatially, enabling si¬ 
multaneous, co-channel transmissions (see e.g. in, El)- While 
it is known that the spatial matched Alter performs optimally 
when the number of base station antennas is large, it is also 
known that there is a signiflcant range of realistic system 
parameters for which the linear MMSE receiver can provide 
spectral efflciencies equivalent to systems with matched-fllter 
receivers with approximately an order of magnitude more 
antennas 0. Hence, it is interesting to study the performance 
of systems with the MMSE receiver and moderately large 
numbers of base-station antennas. 

The uplink of massive MIMO systems with MMSE process¬ 
ing has been analyzed before in 0, 0, and 0. However, 
the spatial distribution of users and base stations was not 
considered in these works. Explicitly modeling spatial user and 
base station distributions can provide insight into the effects of 
tangible system parameters such as the density of base stations 
and users, and the number of antennas, on data rates. 

Modeling the spatial disAibution of mobiles and base sta¬ 
tions for uplink channels is known to be challenging as noted 
in 0. In that work, the complexity associated with modeling 
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the spatial disAibution of mobiles and base stations is handled 
by forming Voronoi cells about the mobiles and disAibuting 
one base station per Voronoi cell. This approach results in 
the Aansmitters in the network being disAibuted according 
to a Poisson Point Process (PPP) for which a number of 
analytical techniques have been developed recently. In other 
works that consider power conAol, which can be used to 
manage interference as well as conserve battery life, Monte 
Carlo simulations are used (e.g. 0, 0). In 0, the authors 
proposed a method for analyzing power control on the uplink 
of cellular networks which is based on approximating the 
interference as originating from sources in concenAic rings 
about a test base station. These approximations that have been 
taken in the literature underscores the difflculty in analyzing 
the uplink of power-conAolled cellular sytems. 

In systems with many antennas at base stations however, 
asymptotic effects can be utilized to address the complexity of 
the uplink, without having to make signiflcant approximations 
on the spatial user and base-station distributions. In 0, we 
considered the uplink of cellular networks with a similar model 
as that used here, except we considered a specific power 
control model in which mobiles attempt to invert the effects 
of the path loss between themselves and their respective base 
stations. Since the locations of the mobiles were assumed to be 
independent in that work, the Aansmit powers of the mobiles 
were also independent, simplifying analysis. In another related 
work ifTOll . we considered the uplink of a cellular network with 
large antenna arrays at Poisson-disAibuted base stations, but 
with a limited number of active mobiles per cell and no power 
control. 

In this work, we consider a cellular network with hexagonal 
cells and multiantenna linear MMSE receivers. We assume a 
general power conAol algorithm where mobiles are assigned 
Aansmit powers as a function of the positions of all mobiles in 
theA respective cells, relative to their base stations. Thus, the 
Aansmit powers of the mobiles in each cell are correlated. Note 
that this correlation increases the complexity of the analysis 
as compared to the system with independent transmit powers 
we considered in 0. Eurther discussion of this is provided 
in Section III-B. We find that when the number of antennas 
at the base station is moderately large and the network is 
interference limited, with a large number of mobiles, the 



Fig. 1. Illustration of a cellular network with at most A' = 10 users 
transmitting with non-zero powers. Mobiles are represented by the dots and 
circles are used to highlight the active mobiles. 

correlation between the transmit powers of mobiles within 
each cell does not signficantly impact the achievable data 
rates. Hence, the transmit power of each mobile does not 
need to depend on other mobiles in its cell. This finding has 
significant practical implications as it simplifies power control 
computation and reduces the need for base-stations to assign 
transmit powers to mobiles in amplitude-reciprocal systems as 
described further in Section III-B. 

II. System Model 

Consider a 2-dimensional network with base stations located 
on a hexagonal grid and where mobiles connect to their closest 
base stations. Let the density of base stations be pc = 1/^c, 
where Ac is the area of each cell. We further assume that there 
is a base station at the origin, with the locations of the base 
stations denoted by Bq, Bi, B 2 , ■ ■ ■ ■ Bq = o shall be referred 
to as the representative base station. Let the cells be denoted 
by Cq,Ci, ■ ■ ■. Suppose that there is a mobile at Xq, which 
we call the representative mobile, located at some point in the 
cell containing the origin. In the remainder of this work, we 
shall analyze the link between the representative transmitter 
and the representative base station which we shall also refer 
to as the representative link. 

Overlaid on the network of base stations is a radius-i? 
circular network centered at the origin as shown in Figure 
[U In this network, n additional mobiles are at independent 
and identically distributed (i.i.d.) locations. These mobiles, if 
active will be co-channel interferers to the representative link. 
Let Xi, X 2 ,■■■ , Xn represent the positions of the mobiles, 
and Vi = \Xi\ be the distance of the i-th mobile from the 
origin. The density of potentially active mobiles pm satisfies 

n=TTpjnR^- ( 1 ) 


We denote the cell in which the i th mobile falls as Ci. 
The average power (averaged over the fast fading) received at 
each antenna of the representative base station from a mobile 
transmitting with power Pi, from a distance r^, is Pir~°‘. The 
path-loss exponent a > 2 is assumed to be a rational number. 
The transmit power of the i-th mobile is a function of the 
location of all mobiles in its cell, relative to the base station 
in that cell. In other words 

Pi = g {{Xj - Center(Ci) : Xj € Ci}) < Pm , (2) 

where Center(Ci) is the location of the base station at the 
center of the cell Ci. Moreover, we require that the mobiles 
either be assigned zero transmit power or a positive transmit 
power that is bounded from below by Pet,. In other words, if 
Pi > 0, then Pi > Pet,. This is a technical requirement that is 
reasonable in practical systems as mobiles can be expected to 
have a minimum power if they are transmitting. 

(|2| describes a general power control algorithm where the 
transmit power of a mobile is a function of the positions of 
all the mobiles in its cell relative to the location of the base 
station. Additionally, we restrict the total number of mobiles 
transmitting in each cell to K by requring that for each i, 

I {j : Xj e Ci and P^ > 0} | < AT. (3) 

This restriction helps to meet quality of service requirements 
and to permit channel estimation using orthogonal pilot se¬ 
quences as has been proposed for massive MIMO systems in 
the literature (e.g. see ||2l). We assume that the representative 
transmitter is always active, i.e. Pq > 0. At a given sampling 
time, the sampled signals at the N antennas of the represen¬ 
tative base station are represented in y G as follows. 

n 

y = To’^goV^To , (4) 

i=l 

where Xi, is the transmitted symbol of the i-th mobile and 
gi G contains i.i.d., zero-mean, unit variance, circularly 

symmetric, complex Gaussian random variables denoted by 
CA/"(0,1). Since we focus on the interference-limited regime 
where the noise power —0, (|4]l does not include noise, which 
means that our results are applicable to networks with a high 
density of mobiles. 

The asymptotic regime considered here is the limit as N, 
n and i? — 00 , such that n/N = c and pm are constant, 
and ([T]i holds. For the remainder of this paper, whenever n, N 
or R ^ 00 , it is assumed that the other two quantities go 
to infinity as well. Note that while our analysis may provide 
insight into the scaling behavior of such systems, we use 
it as a tool to analyze large networks of a fixed size. We 
additionally require lim„_>oo Pr{Pi > 0)n > N which ensures 
that as R ^ 00 , with high probability, there would be a larger 
number of actively transmitting mobiles in the entire network 
than antennas at the representative base station. This ensures 
that the interference covariance matrix R defined below is 



















invertible with probability 1 for all sufficently large n, N, R. 

(5) 


R = Vr-“P,g,g] 


1 

i=i 

We assume that R and go are known at the representative base 
station. Thus, the results here can be interpreted as bounds 
on the performance of real systems which will be subject to 
channel estimation errors which may be significant in certain 
cases when the number of antennas is very large. 

Ill. Asymptotic Spectral Efficiency 
A. Main Results 

Define a normalized version of the SIR at the output of the 
MMSE receiver. 


Pn = N-^^SIR. 
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in probability. Eurther, by applying the bounded-convergence 
theorem mirroring steps in Appendix E in IIT 2 II . we find that 

lim lim |P[log 2 (1+ SIR)]-log 2 (l + SIR)|0. (12) 

c—>-oo N^oo 

Hence, for systems with large antenna arrays at each base 
station but a much larger number of mobiles in the entire 
network than the number of antennas, the spectral efficiency 
7 and its mean are well approximated as follows 
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Using the standard formula for the SIR at the output of the 
MMSE receiver, we can write 


If we now assume that Pq is a random variable, for large 
N,n,R, but with n N, we can approximate the CDE of 
the spectral efficiency as 


(7) 


With these definitions, conditioned tq, Pq, and denoting 
the probability-density function (PDE) of transmit powers by 
fp{p), we have the following theorem. 

Theorem 1: Consider the network model from Section ttn 
As N,n,R —>■ 00 , Pn —>■ /? in probability where /3 is the 
non-negative solution to the following equation 
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( 8 ) 


/ ^ fp{x)x‘-dxdT = 

Jt( ° ) 2 2praT^ 

\ Pm / 

Proof: Given in Appendix lAl 

Moreover, since Pi < Pm, we can bound the second term 
on the LHS of (l 8 ]l as 


fp{x)x <^dxdT < 


(9) 


00 , the integral converges to 
a f 2 'k\ 


( 10 ) 


Additionally, if we assume all mobiles transmit using Gaussian 
codebooks, by the continuous-mapping theorem (e.g. uni 
(3.103)) and removing the normalization of the SIR 


= 0 
( 11 ) 


2n^p^E ypi 

where Fq{(f) is the CDE of g = Porjj"". 
B. Discussion 


Note that the spectral efficiency is primarily a function of 
the transmit power of the representative mobile and E pi 
which means that there is little dependence between the trans 
mit powers of mobiles within a cell which greatly simplifies 
the power control computation. Moreover, in systems which 
use full-duplex operation and where uplink and downlink 
channels have amplitude reciprocity, mobiles can estimate the 
path losses between themselves and their base stations and 
control their transmit powers, without having the knowledge 
of the transmit powers or path-losses of the other mobiles 
in their cell. Thus, mobiles can control their transmit power 
without requiring the base station to communicate information 
to them which results in a reduction in the protocol overhead. 

Einally, it is worth noting that Theorem 1 in this work 
has a very similar form to Theorem 1 in ||3- However, the 
assumptions used here are quite different which requires a 
different proof. In both these works, the main result is a 
consequence of the convergence in probability of the empirical 
distribution function (e.d.f.) of N 2 Pir~°', for i = 1 , 2 , • • • ,n. 
In M, it was assumed that the transmit power of a given 
mobile was solely a function of the distance between itself 
and its given base station. Thus, the transmit powers of the 
mobiles were not coiTelated. As a result, the e.d.f. converges 
as a direct consequence of the law large numbers. In this 
work, the correlation between transmit powers results in a 
more complicated proof of the convergence of the e.d.f. of 
Pir~°‘, for i = 1, 2, • • • , n as given in Appendix lAl 






































Fig. 2. Simulated and asymptotic CDFs of spectral efficiency for N = 8,16, 
32 and 64 antennas, and iC = 10 with con'elated transmit powers. Mobiles 
closer to their base station transmit with power 0.5 and mobiles farther away 
from their base station transmit with unit power as described in Section M\ 

IV. Numerical Simulations and Results 

We conducted Monte Carlo simulations to verify the ac¬ 
curacy of our asymptotic findings. In all cases, we simulated 
1000 realizations of a cellular network according to the system 
model, with path-loss exponent a — A, density of mobiles 
Pm = 1 and the density of base stations pc = 10“^. In 
all cases, the number of mobiles assigned non-zero transmit 
power K is hxed. Thus, we use a high relative density of 
mobiles to base-stations to ensure that with high probability, 
each cell will have K actively transmitting mobiles (i.e. 
mobiles assigned non-zero transmit powers). 

In Figure |2] we show simulated spectral efficiencies of an ad 
hoc power control algorithm which serves mainly to illustrate 
that correlation between transmit powers does not influence 
the spectral efficiency significantly. Mobiles are assigned three 
different transmit powers, Pg = 0.5, Ph = A and zero. First, at 
most K mobiles from each cell are selected to transmit with 
equal probability from the mobiles in each cell. In each cell, 
half of the mobiles closest to the base station are assigned 
power Pi and the other half are assigned Ph- Using order 
statistics and the PDF of tq which is available in ID, we 
can compute the CDF of the spectral efficiency using (fT^ . 
This computation is straightforward but tedious and is not 
included here for the sake of breivity. The resulting CDF is 
illustrated by the solid lines in the figure and the simulated 
values are represented by the markers. The close agreement 
between the theoretical prediction and simulations verifies the 
accuracy of (fl4l i when there is a high degree of correlation 
between transmit powers in a given cell. 

Figures 0 and 0] illustrate systems which use fractional 
power control (see e.g. Q). As described in the previous 
paragraph at most K mobiles transmit with non-zero power 
per cell, and any mobile which is less than unit distance 
from its base station is assigned zero power in order to avoid 
unbounded received signal power at base stations. Of the 
mobiles selected to transmit, a mobile at Xi transmits with 
power \Xi — Center(Ci)|“'^, where e G [0,1]. When e = 0, the 



Fig. 3. Simulated mean spectral efficiency with path-loss inversion power 
control. Markers represent simulated results, solid lines represent the asymp¬ 
totic predictions and the dashed lines are simulated standard deviations of the 
spectral efficiency. The maximum number of actively transmitting mobiles in 
each cell is limited to the values of K indicated in the plot. 

mobiles transmit with equal power and when e = 1, mobiles 
invert the effects of the path losses between themselves and 
their respective base stations. 

Figure [3 illustrates the mean spectral efficiency vs. number 

of antennas when e = 1. Note that since the path-loss is 

inverted in this model, the spectral efficiency of a given link is 

2 _ 

approximately equal to (fOT l. The term E\P“] for this model 
can be obtained by elementary calculus using existing results 
from 191 and is not included here for brevity. Note here that 
the simulated mean spectral efficiencies and the asymptotic 
prediction are extremely close which helps verify the accuracy 
of the asymptotic results. Moreover, the standard deviation 
(and variance) of the spectral efficiency decays with N which 
supports the conclusion that the spectral efficiency approaches 
its asymptote in probability since mean square convergence 
implies convergence in probability. 

Figure|4]illustrates simulated CDFs of the spectral efficiency 
and its asymtotic approximation from (O, with e = 0.5. 
Fq{q) for this model can be computed using straightforward, 
but tedious calculations by starting from the CDF of the 
distance of a random point in a hexagonal cell to the center 
of the hexagon, given in 0- Observe that with increasing N, 
the accuracy of the asymptotic predictions improves, becoming 
virtually indistinguishable from the simulations when the num¬ 
ber of antennas k N — 128. The asymptotic approximation to 
the CDF of the spectral efficiency for e = 0,0.25, 0.5,0.75,1 
and V = 64 is plotted in Figure |5] This figure indicates that 
for low outage probabilities, e = 1, i.e. path-loss inversion is 
the optimal strategy for fractional power control. 

V. Summary and Conclusions 

We analyze the uplink spectral efficiency in spatially dis¬ 
tributed cellular networks with large numbers of base-station 
antennas and power control. Base stations are assumed to use 
the linear MMSE receiver, and transmit powers of mobiles 
within the same cell are correlated. Simple asymptotic ap¬ 
proximations, which are precise in an asymptotic sense, are 
provided for the spectral efficiency of a representative link 












1 



Fig. 4. Simulated and theoretical CDF of the spectral efficiency with frac¬ 
tional power control, where the 2 -th mobile transmits with power | | ^. 

The markers represent simulated values and the solid lines represent the 
asymptotic approximation. The maximum number of actively transmitting 
mobiles per cell is K = 10. 



Fig. 5. CDF of the spectral efficiency based on the asymptotic approximation 
for e = 0,0.25,0.5,0.75,1 and N = 64. Observe that for low outage 
probabilities, e = 1 results in the largest spectral efficiency. 


in such networks. It is also found that when the number of 
base-station antennas is moderately large, and the number of 
mobiles in the entire network is very large, the correlation 
between the transmit powers of mobiles in the same cell does 
not effect the spectral efficiency significantly. As a result, 
simple power control algorithms with no dependency between 
transmit powers of mobiles is sufficient in such networks. This 
finding thus helps reduce the complexity of practical systems 
and simplfies their analysis. In future work, optimal power 
allocation strategies based on this framework shall be explored. 

VI. Acknowledgement 

We thank Dr. Gunnar Peters of Huawei Research for sug¬ 
gesting this work, and the reviewers for helpful comments. 

Appendix 

A. Proof of Main Result 

The main result is proved using Lemma 1 of ini which is 
repeated here for convenience. 


Lemma 1: Consider the quantity 
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In 
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= —si' —S’J'Sl 
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-1 
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(15) 


where s G and S G comprise i.i.d., zero-mean, 

unit-variance entries, n/N = c, and \I/ = diag(^i, ■ • • fn)- 
Suppose that as n, V —>• oo, the e.d.f. of (V’l, i/’ 2 ,' "f’n) con¬ 
verges in probability to H{x). Additionally, assume that there 
exists an no such that Vn > no, Xmin ^ ^ib > 0, 

with probability 1 , where Xmin{A.) denotes the minimum 
eigenvalue of any square matrix A, and Xa is an arbitrary, 
strictly positive number. Then, 7 in probability where 

7 equals the unique non-negative real solution for 7 in 


1 = 7 c 


rdH (t) 

1 + TJ 


Comparing (flSl l to (|2ll and Q, we observe that 


(16) 


j3N=rf^Po7 (17) 


if s = go, the columns of S are gi,g2,'' gn ™cl 

Thus to use Lemma [T] we need to show that the e.d.f. 
of A“/^Pi A“/^P 2 • • • , A“/^P„ converges in 

probability to a limiting function and find its form. Addition¬ 
ally, we need to show that the minimum eigenvalue property 
required in Lemma [T] is satisfied by this model. 

Let Hn{x) denote the empirical distribution function of 
A“/2Pi A“/ 2 P 2 • • • , A“/2p„ r"". The first con¬ 

dition is proved in the following lemma. 

Lemma 2: As n,N,R 00 , Hn{x) converges in proba¬ 
bility to 


_ 2 , 

H{x) = 1 - PLL^Elph 

r -7 



fp{Pj)dPj. 


(18) 


Proof: Given in Appendix\B\ 

The minimum eigenvalue condition is satisfied as shown 
in Appendix O Note that (fT^ has the same form as the 
corresponding equation in la, even though it was assumed 
in that transmit power of the mobiles are independently 
distributed. Following the steps used to prove Theorem 1 in 
||9l and the fact that fp{Pi) — 0 for P^ > Pm , we find that 
Pn —/3 in probability, where /3 is the unique positive real 
solution to dSll. 


B. Proof of Lemma |2] 

Let Hn{x) — E[Hn{x)]. Following steps used to derive 
(38) in lfT3l . we have that for each e > 0 , 
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Let Rc denote the radius of the smallest circle containing a 
given cell in its interior, and Ca denote the union of all cells 
wholly contained in B{0,R), i.e., Ca = \Jk-.CkGB(o,R)^k, 
where B{Y, D) denotes a disk of radius D centered on Y. 
Xi and Xj are in different cells if \Xi — Xj\ > 2Rc. Let the 
event that Xi and Xj are in different cells, and Xi,Xj G Ca 
be denoted by A. Note that as i? —> oo,Pr(^) —>• 1 as 
the probability of Xi and Xj falling in a cell that is not 
wholly contained in the circular network goes to zero, and 
the probability that Xi and Xj are separated by a distance 
less than 2Rc goes to zero as Q{R~^) from lfT4l 

Note that Pi may be dependent on Center (C^) |, where 
Center (Ci) denotes the center of the cell Ci, where recall that 
Ci is the cell containing Xi. Hence, Pi and = \Xi\ are 
correlated in general. However, we note that since Pi is only 
a function of the mobiles in Ci, for R sufficiently large, 

Pr(ri < x\Pi) < Pr(y < x + 2Rc) (20) 

where j/ is a random variable distributed with uniform prob¬ 
ability in Ca- Recall that A is the event that Xi and Xj are 
in different cells which are wholly contained in B{0, R). For 
R large enough such that |Cyi| > Ac we can write an upper 
bound which is based on (l20l l and assuming that Cj is located 

outside B ^0, ^ ■ 



Observing that Pi and Pj are indepeendent since Xi and Xj 
are in different cells, we can take the expectation of the above 
expression with respect to Pi as follows. 



Next, we take the limit as iV, n, i? —5> oo in the manner of 
Theorem 1, and interchange the order of the expectation and 
limit by the bounded convergence theorem. Noting that as 
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Similar to (|2TI) . we can write a lower bound as follows. 
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omparing (l24l i and ( 1211 1 . we note that in the limit as R ^ oo, 
)th the upper and lower bounds will co-incide. Thus, (l23T l 
)lds with equality. Moreover, following a similar set of steps, 
e can find that 
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d Xj and multiplying (|2^ and (1251) . 
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This implies that 
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Substituting (l26T l into ( fT^ . we have that as R ^ oo. 


Pr(|i/„(x)-iJ„(x)| > I) ^0. (27) 


From the definitions of Hn{x) and Hn(x), 
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Substituting (l25T l and recalling that Pr(^) —1, we have 
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The previous expression implies that for each e, there exists 
an i?o, such that MR > Rq 
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When combined with (IZTT i. we have the desired result; 
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C. Minimum Eigenvalue Condition 

Letting T be the set of mobiles whose transmit power is 
non-zero we have: 


1 , iV? ^ j. 


N 


= Pf 


R° 


idT 


lb 


+ iVt-1 ^ P,(r-“ - i?-“)g.gj (32) 

i&T 

— §*§1 +- p~°‘)sisl 


R° 


idT 


ieT 


+ 


N2 




R° 


~ P^b)g^gl 


idT 


+ iVt-1 - Pib){r-^ - i?-“)g,gl, (33) 

idT 


where follows from (39) in 1131 . The smallest eigenvalue 
of the matrix in the brackets in (l3^ was shown to be bounded 
from below by a positive value A^b with probability 1 for 
N sufficiently large, in Lemma 3 of lfT3l . Moreover, since 
Pi > Pib for all i € T, the remaining two matrices in the sum 
on the RHS of (|3^ are non-negative definite. Thus, by the 
Weyl inequality, the smallest eigenvalue of the matrix on the 
LHS of (l3^ is bounded from below by Xa = PaXib, with 
probability 1 for N sufficiently large. 
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